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Abstract. - It is shown that the three-body forces in the I/7/2 shell, for which recently evidence 
was found on the basis of spectroscopic properties of the Ca isotopes and N = 28 isotones, can 
be most naturally explained as an effective interaction due to excluded higher-lying shells, in 
particular the 2^3/2 orbit. 
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■ According to the shell model, the ground state of "^^Ca 
with Z = 20 protons and iV = 20 neutrons has completely 
filled Is, Ip, Id and 2s proton and neutron orbits. Valence 

eutrons from A'^ = 21 on occupy the I/7/2 orbit which is 
closed at the magic number A = 28. Similarly, protons 
ddded to the magic nucleus ^®Ca occupy the proton I/7/2 
orbit. A natural description in the shell model of nuclei in 
t'his mass region is therefore obtained by restricting neu- 
trons and protons to the I/7/2 orbit, an approach which 
"vV^as extensively used in the past [IHl]. From the energy 
levels of nuclei where these neutron and proton orbits are 
filled, it is evident that there must be strong perturba- 
tions from higher configurations. The features predicted 
for relative positions of levels by using two-body effective 
interactions between the protons and between the neu- 
trons, are only roughly obeyed. This is not surprising 
due to the proximity of higher orbits. In ^^Ca, the 3/2~ 
Ifevel, associated with the 2^3/2 orbit, lies only about 2 
MeV above the 7/2~ ground state. Effects of configura- 
tion mixing have been considered many years ago for the 
calcium isotopes [SJ|5] and for the A = 28 isotones [71IH] 
with results that turned out to be in rather good agree- 
ment with experiment. 



Effective interactions are determined by the model 
spaces for which they are intended. If pure (I/7/2)" config- 
urations are adopted, two-body interactions do not yield 
good agreement with experiment. Hence, three-body ef- 
fective interactions have been introduced [9l[T0]. This ap- 
proach may well reproduce level energies in a better way. 
The pure I/7/2 nuclear states, however, will be only an 
approximation to the correct shell-model states. Even if 
effective operators, other than the Hamiltonian may pro- 



vide some help, they cannot provide information on more 
complicated processes like beta decay. Recently, a detailed 
discussion of the effects of three-body interactions on var- 
ious spectroscopic observables was presented by Zelevin- 
sky [11]. 

Recently, three-body interactions were introduced, in 
addition to two-body ones, for dealing with the proton 
and neutron I/7/2 shells [T^. The author considers his 
results as a "manifestation of three-body forces in I/7/2- 
shell nuclei". The good agreement he obtains could be 
evidence of such forces if the I/7/2 shell was a pure one. 
Even within the shell model, whose wave functions are 
not the real ones, the I/7/2 shells are strongly perturbed. 
Such three-body interactions play a role only if the space 
used is a pure I/7/2 shell. It may happen that introducing 
four-body and higher interactions will improve the fit to 
the data [TU] but this will imply a severe loss of predictive 
power. It seems to us that it is more reasonable to consider 
possible interactions with higher configurations. Attempts 
in this direction were made in the past and good agreement 
with experiment has been obtained by using only two- 
body forces. In fact, we aim to show that the effect of 
mixing with nearby configurations leads to a three-body 
effective interaction rather similar to those of ref . [12] . 

The effect of perturbations of the pure I/7/2 shell, due 
to two-body effective interactions, may assume the form 
of additional two-body interactions as well as three-body 
ones. Such effects were considered in the past in atomic 
spectroscopy [T^fTl] . The aim of this note is to derive, by 
adopting a simple approximation, these additional two- 
body and three-body interactions arising from second- 
order perturbations. In pure (7/2)" configurations, all 
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eigenstates have definite seniority for any two-body inter- 
action. There is experimental evidence that the two lowest 
J = 4 levels in ^^Cr are mixtures of states with seniorities 
V = 2 and v = 4. This seniority mixing was obtained in 
ref. |12j by using a three-body interaction. This mixing, 
however, was obtained in refs. [71IH] due to mixing of these 
states with those obtained by raising one I/7/2 proton into 
the 2p3/2 orbit. 

Consider first the j" configuration. The two-body in- 
teraction V is taken to be a perturbation on the single- 
nucleon Hamiltonian. The first-order contribution of V is 
its expectation value in the states considered. In second 
order, the contribution to the energy of the state a J of the 
j" configuration is given by 



E 



{aJ\V\cJ){cJ\V\aJ) 
E{aJ) - E{cJ) '■ 



(1) 



where the summation is over all states cJ which have non- 
vanishing matrix elements with the state aJ . Such states 
are in configurations which differ from that of aJ in single- 
nucleon states of at most two nucleons. These are the 
j"~^j' and j"~'^j'j" configurations. It is well known that 
contributions of the interaction with the latter configu- 
rations yield, in second-order perturbation, effective two- 
body interactions. In the following, the contributions of 
configurations cJ which differ from the aJ configuration 
by the state of one nuclcon will be considered. 

The energies of all states of the j"~^j' configuration cJ 
may differ by small amounts in the difference E{aJ) — 
E(cJ). In such cases it may be a good approximation to 
replace this difference by a constant term, independent of 
the quantum numbers which characterize the state with 
the given J in the excited configuration c. It is then pos- 
sible to carry out the summation over all states cJ. In 
the following we make this approximation and we also as- 
sume that the energy differences in ([T]) are the same for all 
states in the aJ configuration for any number n of (iden- 
tical) nucleons. This is in accordance with the strict rules 
of perturbation theory where these differences are deter- 
mined only by the single-nucleon Hamiltonian. It is a very 
good approximation for electron states in atoms; in nuclei, 
it is not expected to be a very good one. 

Using second quantized operators, explicit expressions 
may be obtained for the resulting two-body and three- 
body operators. The part of the Hamiltonian which con- 
nects the j" and j"~^j' configurations is given by [TS] 



Y.{fj\V\jfJ)A\jfJM)A{jjJM), (2) 

JM 

where A{jfJM) = {A^jj'JM)y and 



From the formula ([T]), applied to the intermediate config- 
urations cJ of the type jj' J and assuming equal denomi- 
nators, the following sums are obtained: 

E {fJi\Vbj'Ji)A^ (jj Ji Afi) AO'/ JiMi) 

Ml Ah 

X {fj2\V\jj'J2)A^ijfJ2M2)A{jjJ2M2). (4) 

With the help of the anti-commutation relations between 
creation and annihilation operators as well as some tensor 
algebra identities, this expression can be evaluated. The 
result is the sum of a two-body interaction 

E \{fJ\V\jfJ)\'A\jjJM)AijjJM), (5) 

JM 



and a three-body one 
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J1J2JM 



{fjl\V\jj'Jl){fj2\V\jj'J2){ ] jI ^j } 



X [AHjjJi) X a]](^^^) [A(.7j J2) X a,]^-'''\ (6) 

Instead of using these expressions, it is simpler to reduce 
the matrix elements in the j" configuration to those in the 
p and j^ configurations, where the calculation of matrix 
elements is straightforward. These matrix elements define 
uniquely the operators. To calculate matrix elements be- 
tween states aJ and hJ in the j" configuration one may 
use the expansion |15| 

(7) 
where [j" ^ {aoJojj J\} j^ aJ] are coefficients of fractional 
parentage (c.f.p.) which express the anti-symmetric n- 
particle wave function "^(j^aJ) in terms of an anti- 
symmetric (n—l)-particle wave function 'I'(j"~^ao Jo) cou- 
pled with that of the n*^ particle j„ to total angular mo- 
mentum J. Since this wave function, as well as the one 
of the bJ state, is fully antisymmetric, it is possible to 
replace the matrix element of any two-body interaction 
^ = Sife ^fc by V12 multiplied by n{n — l)/2, the number 
of such terms. Thus, we obtain 



irajinrbj) 

n{n — 1) 

2 
n{n — 1) 



{raJ\Vi2\rbJ) 

{faoMV^2\r-%.h) 



(3) 



2 ^^ 

aoboJo 

Mr'\aoJo)jj\}raj][r-\boJo)jj\}rbj] 
= ^ E {r~'^oJo\v\r-'boJo) 

aoboJo 

x[r-\aoJo)jJ\}raJ][r-\boJo)jJ\}rbJ].{8) 

This procedure may be further applied to the matrix ele- 
ments in the j"~^ configuration until the j^ configuration 
is reached. 
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Similarly, it is possible to apply this procedure to any To evaluate the first sum it is convenient to carry 

three-body interaction V = X^iifc ^'jk 8'iid to replace it out a change of coupling transformation on the state 

by n{n - l){n - 2)Vi23/6. Using c.f.p., it is possible to iJshWfiJ) = -\hh{Ji).fiJ), 

calculate matrix elements in the j" configuration in terms .^^ . ., „ , 

of matrix elements in the j"~^ configuration. Thus, we \j2j3[Ji)JiJ} = 2^^^ ' J1J2 

obtain -^^ 

{raJ\V\rbJ) ^ { J ' jI j- }b2jI(J2)j3J),(13) 

_ n(n-l)(n-2) y , .„, , 

- Q \-? a-i\yi23\J OJ) ^YiBYe J, = ^f^~T\. The integration over J3 can be 

nin — l)(n — 2) v-^ - carried out and yields the non-vanishing terms in which 

== ^ 2^ 0" aidJQ\y\2'i\r bo Jo) the equality J2 = Jo must hold, as follows 

aaboJo 

X [.r-\aoJo)jj\}.raj][r-\boJo)jj\}rbj] vsi^iy^^' ^ JoMfWjJiyfaj] 



J2 {.r''<^oJo\v\r-%Jo) 



Jo 



" ^-""^o' x{ ^; -^0 •( }{fJo\V\j2j[Jo). (14) 

Mr-\aoJo)jJ\}raJ][r-\boJo)jJ\}rbJ]. (9) ^ J Ji J J 

Similarly, the second summation in (J12p yields the follow- 

This procedure may be continued until the j configura- • i, . 

tion is reached. 

Let us now turn to the calculation of the effective inter- —v^f—i U+i' \^ / / \j^( j )n jlXj^n 71 

action in the j^ and j^ configurations, due to additional , 

shells j' that are not explicitly taken into account. The j j 

two-body operators may be obtained from perturbation x -^ '., ({J^Jol^lJiJ-zJo}- (15) 
theory in the two-nucleon configuration. The contribu- 
tion in second-order perturbation to the f configuration ^j^g ^^^ expressions may be combined yielding the result 
in a state with spin Jg can be simply calculated from ([T]) 
and is given by _V6(-1)^+^' ^ JoMf{Jo)jJ\}faJ] 

2-.E{pjo)-E{jj'Joy ^^^^ x{J, 2 ,Yf-Mv\n'.h), (16) 

where the matrix elements of the interaction V are be- where again the states {pJo\ and |jj'Jo) are anti- 

tween anti-symmetric and normalized states (jVo| and symmetric and normalized. To obtain the second-order 

\ij'Jo). contribution ([T|) of these perturbations, this expression 

To obtain the matrix elements of the three-nucleon con- should be multiplied by a similar one where the summa- 

tribution, the f configuration should be considered. Re- tion is over all values of J^. The second-order perturbation 

suits of perturbations of states in the j^ configuration may contribution in the f configuration is therefore 

be calculated as follows. Anti-symmetric and normalized .^^ .^^ ^ ^ 

cJ states may be expressed as " Z^( ^ ^ -' Z^ '^oJq 

j'.h .laJ'a 

nfwfj) = ^^{nfwfsJ)-nf{Ji)j[j) x\j\jo)jj\}faj]{^, ;;•; ;^ } 

-*(.?'(^i)j2^))- (11) r -i V J 1 



The matrix element of y = V12 + F13 + V23 between the 

fully anti-symmetric states \aJ) = \faJ) and \cJ) = ^, (J Jo\y\JJ' Jg){JJ' Jo\y\J Jp 
\fiJi)j'J) is equal to the matrix element of W12 and E{paJ) - E{f{Ji)j' J) 

hence is given by 



As explained above, we make the approximation that 
/rv^r -2/ T \ ■ rn -3 ti; -2/ t \ ■ tit'a \ -2 r t \ ■/ t\ the energy denominators in (flj are the same for all states 
^_^vj \ Kijj u.y jw V yjjjo i iz|j v i/Ji / ^f ^j^g ^n ^^^ ^n ly configurations. This means that 

" for the y"^ configurations we assume that the energies 

-VsY,[fWjJ\}.faJ]{fiJo)J3J\Vi2\f{Ji)]2J)- E{faJ) = E{f) and E{f{Ji)j'J) = E{fj'), are in- 

•■'o dependent of a, Ji and J, and given be 3ej and 2ej + ej' , 

(12) respectively, where ej and e^/ are single-particle energies. 
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Table 1: The T = 3/2 three-body interaction matrix elements in the I/7/2 shell obtained by Volya and the effective three-body 
interaction matrix elements derived with the formula (|21|l for various p3f7 interactions and for the gxpfla interaction. Energies 
are in keV. 







Z = 20 






TV = 28 






J 


Volva jl2, 


p3f7 |5J 


p3f7 16J 


Volva 1121 


p3f7 |7j 


p3f7 18J 


gxpfla 


3/2 


-559 (273) 


-412 


-173 


-128 (88) 


-53 


-115 


-41 


5/2 


2 (185) 


-207 


-104 


-18 (70) 


-54 


-57 


-75 


7/2 


53 (70) 


138 


76 


55 (28) 


42 


37 


45 


9/2 


272 (98) 


392 


157 


122 (41) 


36 


110 


105 


11/2 


51 (130) 


29 


46 


102 (43) 


72 


6 


79 


15/2 


-24 (73) 


77 


28 


-53 (29) 


5 


23 


-1 



In this case, the summation over Ji can be carried out 
directly due to the identities of Racah coefficients. The 
sum over even values of Ji can be expressed as a linear 
combination of sums in which the summation is over all 
values of Ji , yielding 



J I even 

.E(-i)^n2.t + i){^ t]]{]^l]] 



Sjo 



2Jn 



f 1 + I J f J'^ J 



(18) 



Substituting this result into the expression above we find 
3E E J^4[f{J^)jJ\}faJ][f{4)]J\}3'aJ] 

J' JaJ'a 

5.1„.1L ( i J Jr. 



"JoJq , r 3 J Jo \ 
2Jo + l I J f Jo i 

{fJo\V\jj'Jo){jfJ[>\V\fJl^) 
E{P) - E{pf) 



(19) 



To obtain the matrix elements of the genuine (but effec- 
tive) three-body interaction, the contribution of the two- 
body interaction should be subtracted. The latter is given 
by 



S[fiJo)jJ\}faJrY. 



\{fJo\Vi2\jfJo)\' 
^E{pJo)^E{jj'Jo)- 



(20) 



Again we assume that the energies E{j^Jq) = E{p) = 2ej 
and E{jj'Jf)) = E{jj') = ej + eji , are independent of Jo- 
in that case, subtracting this contribution from (J19p . we 
find the genuine three-body matrix element to be equal to 



3E E JoJo[f{Jo)jJ\}faJ][j\j[,)jJ\}faJ] 

j' JoJ'o 

^ r J J Jo \ {3^Jo\v\n'Jo){jfJo\v\fJ'o) 

I J f Jl, / E{3^) - E{3^r) 



(21) 



It should be pointed out that the second-order pertur- 
bation contribution is contained in both two-body and 
three-body interactions. The contribution of second or- 
der in perturbation theory is always attractive (negative). 
This need not be the case for the three-body interaction 
which may lead to some strange results. For example, if 
j = 7/2 and j' = 1/2, there is no contribution of this 
perturbation to the j^ state with J = 15/2. The highest 
value of Ji is Ji = 6 which cannot couple with j' = 1/2 
to yield J = 15/2. Still, there is a non- vanishing contribu- 
tion to the matrix clement with J = 15/2 and Jq = 6 of 
the two-body interaction. This unphysical contribution is 
exactly cancelled by the corresponding contribution from 
the three-body interaction. Such cancellation must take 
place in all 7/2" configurations. 

In the calculation of matrix elements in the j^ configura- 
tion, the origin of the unphysical terms is clear. The sum 
of products of two 6j-symbols in (|18p is always l/(2Jo + l) 
even if all 6j-symbols with Ji even vanish. In such a case, 
the result is due to symbols with odd Ji values. These 
symbols contribute to a non- vanishing symbol in which Ji 
no longer appears. If in ((T5)) Jq ^ Jg, there is no con- 
tribution of terms with odd Ji values to the three-body 
interaction. If, however, Jq = Jq, the odd Ji terms con- 
tribute with opposite signs the same amount to the two- 
body interactions and to the three-body ones. Thus, the 
two-body (PO)) and three-body (PT|) expressions calculated 
above, which are genuine two-body and three-body inter- 
actions, may contain unphysical terms. The contributions 
of two-body interactions are usually absorbed into the ef- 
fective two-body interaction. If the interest is in three- 
body interactions whose contributions cannot be mim- 
icked by two-body terms, the expressions obtained above 
should be used. 

The above formalism can be applied to the I/7/2 shell. 
Several interactions are available that include effects from 
the 2p3/2 shell. Two of them were derived on the basis 
of spectroscopic properties of the calcium isotopes [5l|6] 
and two more from those of the A^ = 28 isotones [7l[8] . In 
particular, the authors of these references give numerical 
values for the matrix elements involving the 2^3/2 shell 
which enter the expression ([?!]) and for the difference in 



single-particle energies, £2^3 



=1/7/: 



, . We refer to these in- 
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teractions as p3f7, followed by the relevant reference. For 
comparison, we also include results obtained with a more 
recent interaction which considers the 2plf shell and re- 
sults from an empirical fit (starting from a microscopically 
derived set of matrix elements) to reproduce a large body 
of energy data in the 2pl/-shell nuclei. This is the gxpfla 
interaction (16j which adopts the single-particles energies 
eif^^^ = -8.6240, e2p^^^ = -5.6793, ezpi/o = -4.1370 and 



^1/5/ 



-1.3829 MeV. 



1/2 



We have used these interactions, without any modifi- 
cation, to calculate the effective three-body interaction 
which is induced when the model space is restricted to 
1/7/2- The results are shown in Table [TJ For j = 7/2 
a three-nucleon state is completely specified by its total 
angular momentum J and there is no need for the ad- 
ditional label a. There are six three-nucleon states with 
J = 3/2, 5/2, 7/2, 9/2, 11/2 and 15/2, and each of these 
states defines a component of the three-body interaction. 
According to the expression (|21|) , the total effective three- 
body interaction results from additive contributions of the 
different shells j' which can be 2pi/2i 2^3/2 or 1/5/2- The 
three-body matrix elements of Volya [12] result from sep- 
arate fits to the calcium isotopes and to the A^ = 28 iso- 
tones, and they should thus be compared with correspond- 
ing matrix elements obtained with the p3f7 interactions 
from rcfs. [Hlini a-nd [71[5], respectively. The numbers in 
parentheses in Table [T] are also taken from Volya and cor- 
respond to the variances of the parameters in the fits. 

There are substantial variations in the calculated ef- 
fective three-body matrix elements. In particular, those 
derived from the Z ~ 20 p3f7 interactions are generally 
larger than those obtained with the TV = 28 p3f7 inter- 
actions. This results from a combination of a larger dif- 
ference £2pa/2 ~ £1/7/2 a^^ smaller two-body matrix ele- 
ments ((l/7/2)^-'^o|V"|l/7/22p3/2Jo) in the latter interac- 
tions. Note also that the effective three-body matrix el- 
ements derived from gxpfla are closer to those obtained 
with the iV = 28 p3f7 interactions. 

There are certainly differences between these calcula- 
tions and the results of Volya. However, in spite of these 
differences and uncertainties, it is clear that the matrix el- 
ements are correlated: with one exception, attractive (re- 
pulsive) matrix elements in the analysis of Volya turn out 
to be attractive (repulsive) in our analysis. The exception 
concerns the J = 15/2 interaction which is attractive in 
Volya's analysis while it is repulsive for the p3f7 interac- 
tions and essentially zero in gxpfla. 

It is believed that there are important ab initio three- 
body interactions and also three-body interactions due 
to short-range correlations between nuclcons. The lat- 
ter arise from admixtures of highly excited configura- 
tions. So far, no evidence was found for the effects of 
three-body interactions on states of valence nucleons. In 
cases where rather pure shell-model configurations were 
observed, states and energies were well determined by ef- 
fective two-body interactions. Clearly, the I/7/2 shell for 



protons and for neutrons is not in this category since there 
are low-lying configurations whose states may well mix 
with the (1/7/2)" states. Insisting on using pure I/7/2 con- 
figurations, better agreement with experiment is obtained 
by incorporating effective three-body interactions. We do 
not question the existence of such three-body interactions 
but we surmise that most likely they arise from renormal- 
ization effects due to admixtures of rather low-lying con- 
figurations. Our conclusion is that claims of evidence for 
the existence of genuine three-nucleon interaction should 
be treated with some skepticism if they are based on cal- 
culations in a shell-model space which seems to be too 
restricted. 
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